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A Bethe-Salpeter treatment of Cooper pairs (CPs) based on an ideal Fermi gas (IFG) "sea" 
yields the familiar negative-energy, two-particle bound-state if two-hole CPs are ignored, but is 
meaningless otherwise as it gives purely-imaginary energies. However, when based on the BCS 
ground state, legitimate two-particle "moving" CPs emerge but as positive-energy, finite-lifetime 
resonances for nonzero center-of-mass momentum, with a linear dispersion leading term. Bose- 
Einstein condensation of such pairs may thus occur in exactly two dimensions as it cannot with 
quadratic dispersion. 
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Shortly after the publication of the BCS theory 
of superconductivity charged Cooper pairs (CPs) ob- 
served in magnetic flux quantization experiments with 
3D conventional QQ, and much later with quasi-2D 
cuprate L S\ superconductors, suggested CPs as an indis- 
pensable ingredient. Although BCS theory admits the 
presence of Cooper "correlations," several boson-fermion 
(BF) models |f|-|l£| with real, bosonic CPs have been 
introduced after the pioneering work of Refs. [lfij-[lfij. 
However, with one exception [fj-^l]], all such models ne- 
glect the effect of two-hole (2h) CPs treated on an equal 
footing with two-particle (2p) CPs — as Green's functions 
|20| can naturally guarantee. 

The BCS condensate consists of equal numbers of 2p 
and 2h Cooper correlations; this is evident from the per- 
fect symmetry about fi, the electron chemical potential, 
of the well-known Bogoliubov [2lJ v 2 (e) and u 2 (e) coef- 
ficients [see just below (0 later on] , where e is the elec- 
tron energy. Some motivation for this Letter comes from 
the unique but unexplained role played by hole charge 
carriers in the normal state of superconductors in gen- 
eral [23, as well as from the ability of the "complete (in 
that both 2h- and 2p-CPs are allowed in varying propor- 
tions) BF model" of Refs. to "unify" both BCS 
and Bose-Einstein condensation (BEC) theories as spe- 
cial cases. Substantially higher T c 's than BCS theory are 
then predicted without abandoning electron-phonon dy- 
namics. Compelling evidence for a significant presence of 
this dynamics in high-T c cuprate superconductors from 
angle-resolved photoemission spectroscopy data has re- 
cently been reported 

In this Letter the Bethe-Salpeter (BS) many-body 
equation (in the ladder approximation) treating both 
2p and 2h pairs on an equal footing is used to show 
that, while the ordinary CP problem [based on an ideal 
Fermi gas (IFG) ground state (the usual "Fermi sea")] 



does not possess stable energy solutions: i) CPs based 
not on the IFG-sea but on the BCS ground state sur- 
vive as positive energy resonances; ii) their dispersion re- 
lation in leading order in the total (or center-of-mass) 
momentum (CMM) HH = 7i(k, + k^) is linear rather 
than the quadratic h 2 K 2 /4m of a composite boson (e.g., 
a deuteron) of mass 2m moving not in the Fermi sea 
but in vacuum; and iii) this latter "moving CP" solu- 
tion, though often confused with it, is physically distinct 
from another more common solution sometimes called 
the Anderson-Bogoliubov-Higgs (ABH) 0, (0 p. 44), 
|25(-|28| collective excitation. The ABH mode is also lin- 
ear in leading order and goes over into the IFG ordi- 
nary sound mode in zero coupling. A new feature emerg- 
ing from our present 2D results, compared with a prior 
3D study outlined in Ref. H^], is the imaginary energy 
term leading to finite-lifetime CPs. We focus here on 
2D because of its interest [? ] [? ] for quasi-2D cuprate 
superconductors. In general, our results will be crucial 
for Bose-Einstein condensation (BEC) scenarios employ- 
ing BF models of superconductivity, not only in exa ctly 
2D as with the Berezinskii-Kosterlitz-Thouless |30|[3~uj 
transition, but also down to (1 + e)D which characterize 
the quasi- ID org ano-metallic (Bechgaard salt) supercon- 
ductors |s^-|34j. Striking experimental confirmation of 
how superconductivity is "extinguished" as dimension- 
ality d is diminished towards unity has been reported 
by Tinkham and co-workers pjjj |36l] . They measured 
resistance vs. temperature curves in superconducting 
nanowires consisting of carbon nanotubes sputtered with 
amorphous MojgGe2i and of widths from 22 to 10 nm, 
showing how T c vanishes for the thinnest widths. Our 
results also apply, albeit with a different interaction, to 
neutral-atom superfluidity as in liquid 3 He [13] as well as 
to ultracold trapped alkali Fermi gases such as 6 Li 
and 40 K [^3 since pairing is believed to occur there also. 



2 



For bosons with excitation energy ek = C S K S +o(K s ) 
(for small CMM K) BEC occurs in a box of length L 
if and only if d > s, since T c = for all d < s. The 
commonest example is s — 2 as in the textbook case of 
ordinary bosons with sk — h 2 K 2 /2m exactly, giving the 
familiar result that BEC is not allowed for d < 2. The 
general result for any s is seen as follows. The total boson 
number is 

N = N (T) + [exp/3(e K - Ms) - 1] 1 

K^O 

with (3 = ksT. Since Nq(T c ) ~ while the boson chemi- 
cal potential fis also vanishes at T = T c , in the thermo- 
dynamic limit the boson number density becomes 

I* oo 

N/L d ~ A d / dKK d - l [e^p c {C s K s + •••)- l]" 1 

where Ad is a finite coefficient. Thus 

N/L d ~ A d (fc s T c /C s ) / di^-*- 1 + / • • • , 

where if max is small and can be picked arbitrarily so long 
as the integral J K • • • is finite, as is N/L d . However, 
if d = s the first integral gives InK \^ rnM = — oo; and 
if d < s it gives l/(d - s)K s ~ d |f m ~= -oo. Hence, T c 
must vanish if and only if d < s, but is otherwise finite. 
This conclusion hinges only on the leading term of the 
boson dispersion relation ek- The case s — 1 emerges in 
the CP problem to be discussed now. 

In dealing with the many-electron system we assume 
a BCS-like electron-phonon model s-wave inter-electron 
interaction, whose double Fourier transform j/(|ki — k'J) 
is just 

u(k 1 ,k[) = -(k F /k' 1 )V (1) 

if kp — kD < ki < kp + kp,, and = otherwise. Here 
V > 0, hkp = mvp the Fermi momentum, m the effective 
electron mass, v F the Fermi velocity, and fco = ujd/vf 
with u>d the Debye frequency. The usual condition 
hiou -C Ep then implies that kn/kF = huJo/^Ep <C 1. 

The BS wavefunction equation |2J| in the ladder ap- 
proximation with both particles and holes for the original 
IFG-based CP problem using (JIJ leads to an equation for 
the wavefunction ip-^ in momentum space for CPs with 
zero CMM K = k : + k 2 = that is 

(2&-£b)tk = ^53'V^-VX)"^'- ( 2 ) 

k' k' 

Here ^fc = H 2 k 2 /2m — Ep, £q is the eigenvalue energy 
and k = ^(ki — k2) is the relative wavevector of a pair. 
The single prime over the first (2p-CP) summation term 
denotes the restriction < < hwp, while the double 
prime in the last (2h-CP) term means —frbjp, < < 0. 



Without this latter term we have Cooper's Schrodinger- 
like equation 0] for 2p-CPs whose implicit solution is 

clearly ijjk = (2^ — £o)~ 1 Vj2k> Vv- Since the sum- 
mation term is constant, performing that summation 
on both sides allows canceling the ^k-dependent terms, 

leaving the eigenvalue equation ^ k (2£fc — fo) _1 = 1/V 
with the familiar solution £q = — 2hwp)/ {e 2 / x — 1) (ex- 
act in 2D, and to a very good approximation otherwise 
if huj D < Ep) where A = VN(E F ) with N(E F ) the elec- 
tronic density of states (DOS) for one spin. This corre- 
sponds to a negative-energy, stationary-state bound pair. 
For K ^ the CP eigenvalue equation becomes 

Y,\^k + h 2 K 2 /2m-£ K )- 1 =1/V. (3) 

k 

Note that a CP state of energy £k is characterized only 
by a definite K but not definite k, in contrast to a "BCS 
pair" defined with fixed K and k (or equivalently def- 
inite ki and k2). Without the first summation term in 
(J2J the same result in £q for 2p-CPs follows for 2h-CPs 
(apart from a sign change). However, using similar tech- 
niques to solve the complete equation (0) — which can- 
not be derived from an ordinary (non-BS) Schrodinger- 
like equation in spite of its simple appearance — gives the 
purely- imaginary £q — zti2hcjp,/Ve 2 / x — 1, thus imply- 
ing an obvious instability. This was reported in Refs. 
[24( p. 44 and ^(J wn0 did not stress the pure 2p and 2h 
cases just discussed. Clearly then, the original CP pic- 
ture is meaningless if particle- and hole-pairs are treated 
on an equal footing as consistency demands. This is per- 
haps the prime motivation for seeking a new unperturbed 
Hamiltonian about which to, e.g., do perturbation the- 
ory. 

A BS treatment not about the IFG sea but about the 
BCS ground state vindicates the CP concept. This sub- 
stitution might seem an artificial mathematical construct 
but its experimental support lies precisely in Refs. 
and its physical justification lies in recovering two ex- 
pected results: the ABH sound mode as well as finite- 
lifetime effects in CPs. In either 3D 29] or 2D the BS 
equation yields two distinct solutions: the usual ABH 
sound solution and a highly nontrivial "moving CP" so- 
lution. The BS formalism gives rise to a set of three 
coupled equations, one for each (2p, 2h and ph) channel 
wavefunction for any spin-independent interaction such 
as Q. However, the ph channel decouples, leaving only 
two coupled wavefunction equations for the ABH solu- 
tion. The equations involved are too lengthy, and will 
be derived in detail elsewhere. The ABH collective exci- 
tation mode energy £r is found to be determined by an 
equation that for K = gives £ = (Ref. p. 39) 
and reduces to J ° d^/^/^ 2 + A 2 = 1/A, the familiar 
BCS T = gap equation for interaction whose solu- 
tion is A = hwp>l sinh(l/A). Taylor-expanding £k about 
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K = and small A gives 

flVF 



£ 



K 



K + Q{K 2 ) 



V2 



(4) 



Note that the leading term is just the ordinary sound 
mode in an IFG whose sound speed c = v F /vd in d 
dimensions which also follows trivially from the zero- 
temperature IFG pressure P = n 2 [d(E/N)/dn] = 
2nE F / (d + 2) on applying the familiar thermodynamic 
relation dP/dn — mc 2 . Here E — dE F /(d + 2) 
is the IFG ground-state energy while n = N/L d — 
kp/d2 d ~ 2 TT d / 2 T(d/2) the fermion-number density. 

The second solution in the BCS-ground-state-based BS 
treatment is the moving CP solution for the pair energy 
£k which in 2D is contained in the equation 

—XHvp / dk difu K/2+ kV K /2-kX 

Zn Jk F -k D JO 

x{w K /2-k«K/2+k - UK/2+k«K/2-k}x 

= 1, (5) 



-E-K/2+k + E K / 2 -k 



^K/2-k) 2 




where tp is the angle between K and k; A = VN(E F ) 
as before with N(E F ) = m/2-Kti 2 the constant 2D DOS 
and V the interaction strength defined in QJ; Ey = 
A 2 with A the fermionic gap; while u\ = |(1 + 
k) and v\ = 1 — u\ are the Bogoliubov functions 
[21|. In addition to th e pp and hh wavefunctions (de- 
picted grafically in Ref. [23 Fig. 2), diagrams associated 
with the ph channel give zero contribution at T = 0. A 
third equation for the ph wavefunction describes the ph 
bound state but turns out to depend only on the pp and 
hh wavefunctions. Taylor-expanding £k in powers of K 
around K = 0, and introducing a possible damping factor 
by adding an imaginary term — iT k in the denominator, 
yields to order K 2 for small A 



±£ 



K 



2A H hv F K H 

2tt * 9 k D 



F e 1 ' x K 2 



A 



-hvpK - 



1 Hv f ^ /Xk2 



12 ki 



Q{K 3 



(6) 



where the upper and lower sign refers to 2p- and 2h-CPs, 
respectively. A linear dispersion in leading order again 
appears, but now associated with the bosonic moving 
CP. The positive-energy 2p-CP resonance has a lifetime 
t k = h/2T K = fr/2 [(\/ir)hv F K + (hv F /12k D )e^ x K 2 ] 
diverging only at K = 0, and falling to zero as K in- 
creases. Thus, "faster" moving CPs are shorter-lived 
and eventually break up, while "non-moving" ones are 
stationary states. The linear term (\/2n)hv F K con- 
trasts sharply with the coupling-independent leading- 
term in £ K = £ - {2/Tr)hv F K + 0{K 2 ) (or 1/2 in 
3D instead of 2/tt) that follows from the origi- 
nal CP problem ||3J| neglecting holes — for either inter- 
action Q E| 0T an attractive delta inter-fermion po- 
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FIG. 1: Exact "moving Cooper pair" energy £k (in units of 
Ep) from JSJ (full curves), compared with its linear leading 
term (short-dashed lines) and its linear plus quadratic expan- 
sion (long-dashed curves) both from ©, vs CMM wavenum- 
ber K (in units of kp), for interaction (1) parameters A = 
2 (lower set of curves) and | (upper set of curves), and 
hu)p/Ep = 0.05. For reference, the leading linear term (4) 
of the ABH sound mode is also plotted (lower thick dashed 
line) . 



tential |43] ; 44| (imagined regularized |4Ej to have a sin- 
gle bound state whose binding energy serves as the cou- 
pling parameter). In the latter simple example, more- 
over, it is manifestly clear in 2D that the quadratic 
H 2 K 2 /4:Tn stands alone as the leading term for any cou- 



pling only when E F 



is strictly zero, i.e., in the 



2, giving 
0.004 and 



absence of the Fermi sea. Fig. 1 graphs the exact mov- 
ing CP (mCP) energy extracted from JSJ, along with 
its leading linear-dispersion term and this plus the next 
(quadratic) term from (jSJ. The interaction parameter 
values used in were Hu>d/E f = 0.05 (a typical value 
for cuprates) and the two values A = \ and 
for £ /E F = 2A/E F = 2?vjQ D /E f sinh(l/A) ~ 
0.028, respectively (marked as dots in the figure). Re- 
markably enough, the linear approximation (thin short- 
dashed lines in figure) is better over a wider range of 
K/k F values for weaker coupling in spite of a larger and 
larger partial contribution from the quadratic term in 
©; this peculiarit y al so emerged from the ordinary CP 
treatment of Ref. [42j and might suggest the expansion 
in powers of K to be an asmyptotic series that should 
be truncated after the linear term. For reference we also 
plot the linear term hv F K/y/2 of the sound solution (4). 

We cannot presently address such matters as the na- 
ture of the normal state, the pseudogaps observed in un- 
derdoped cuprates, etc., but efforts in these directions 
are in progress. 

Like Cooper's [see Eq. ©], our BS CPs are charac- 
terized by a definite K and not also by definite k as the 
pairs discussed by BCS @|. Hence, the objection does not 
apply that CPs are not bosons because BCS pairs with 
definite K and k (or equivalently definite ki and k2) have 
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creation/annihilation operators that do not obey Bose 
commutation relations [Ref. 1], Eqs. (2.11) to (2.13)]. 
In fact, either or (J5J shows that a given "ordinary" 
or BS CP state labeled by either K or £^ can accommo- 
date (in the thermodynamic limit) an indefinitely many 
possible BCS pairs with different k's. This implies BE 
statistics for either ordinary or BS CPs as each energy 
state has no occupation limit. 

To conclude, hole pairs treated on a par with electron 
pairs play a vital role in determining the precise nature 
of CPs even at zero temperature, only when based not 
on the usual ideal-Fermi-gas (IFG) "sea" but on the BCS 
ground state. Treatment them with the Bethe-Salpeter 
equation gives purely-imaginary-energy CPs when based 
on the IFG, and positive-energy resonant-state CPs with 
a finite lifetime for nonzero CMM when based on the 
BCS ground state — instead of the more familiar negative- 
energy stationary states of the original IFG-based CP 
problem that neglects holes, as sketched just below J2J. 
The BS "moving-CP" dispersion relation is gapped by 
twice the BCS energy gap, followed by a linear leading 
term in the CMM expansion about K = 0. This linearity 
is distinct from the better-known one associated with the 
sound or ABH collective excitation mode whose energy 
vanishes at K = 0. Thus, boson-fermion models assum- 
ing this CP linearity for the boson component instead of 
the quadratic h 2 K 2 /Am can give BEC for all d > 1, in- 
cluding exactly 2D, and thus in principle apply not only 
to quasi-2D cuprate but also to quasi-lD organo-mctallic 
superconductors. 
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